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This paper is a generalization of the models considered in [J. Stat. Phys. 128,1365 
(2007)]. Using an analogy with free fermions, we compute exactly the large deviation 



\ function (ldf) of the energy injected up to time t in a one-dimensional dissipative 

q system of classical spins, where a drift is allowed. The dynamics are T = asym- 

o : 

metric Glauber dynamics driven out of rest by an injection mechanism, namely a 

Poissonian flipping of one spin. The drift induces anisotropy in the system, making 
OS ' 

' the model more comparable to experimental systems with dissipative structures. We 

\o ' 

1/"} ■ discuss the physical content of the results, specifically the influence of the rate of the 

Poisson injection process and the magnitude of the drift on the properties of the ldf. 
OO ! 

, We also compare the results of this spin model to simple phenomenological models 



of energy injection (Poisson or Bernoulli processes of domain wall injection). We 
show that many qualitative results of the spin model can be understood within this 
simplified framework. 
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I. INTRODUCTION 

Dissipative systems are generically systems for which a few relevant degrees of freedom can 
be singled out and obey closed dynamical equations: typically a fluid, where the velocity field 
obeys the Navier-Stokes equation, belongs to this category. Another well-known example is 
given by granular materials, where the identification of relevant variables (collisions) is even 
more evident. The lack of completeness, caused by the selection of some degrees of freedom, 
gives however these systems a nonconservative character, as energy flows continuously from 
relevant degrees of freedom (kinetic energy) to irrelevant ones (thermal agitation). As a 
result, the dissipative systems are by nature very different from the systems usually suitable 
for the use of classical statistical physics, where the conservation of energy is an unavoidable 
assumption. In particular, the whole set of tools devised by statistical physics can be of 
questionable use, even in situations where a statistical approach seems natural: it is very 
tempting to interpret turbulent systems, or a vibrated granular matter, in terms of effective 
temperature, correlations, Boltzmann factor, etc. . . but the soundness of such an approach 
is often questionable. 

Quite recently, the interest of physicists has been drawn to the injection properties of 
dissipative systems for several reasons. First, it was easily measurable experimentally, and 
the measurements showed that, contrarily to what was usually expected, the injected power 
fluctuates a lot, is not Gaussian, and does not obey the usual simple scaling arguments 

aa. 

Moreover, the injection is by nature very important in dissipative systems, since it is required 
to draw the system out of rest; thus, it is natural to study specifically this observable, which 
is at the same time responsible for the existence of the stationary state, and is strongly 



affected by it 



4|. Finally, some theoretical works on the so-called "Fluctuations Theorems" 



had suggested a possible symmetry relation in the distribution of the fluctuations of the 
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injected power, a suggestion vigourously debated since the works of 

In studying the fluctuations of global (macroscopic) variables of a disordered (turbulent) 
dissipative system, one faces soon a crucial problem: contrary to the statistical physics of 
conservative systems, no global theory is at hand here to predict the level of fluctuations, 
the physical meaning of their magnitude, the skewness of the distributions, etc. . . All these 
features are intimately connected to the statistical stationary turbulent state, but in a way 
nowadays beyond our knowledge. A way to make progress towards a better understanding 
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of these issues is to consider toy-models of dissipative systems where some features of real 
systems are reproduced, and analyze the structuration of the stationary states. If one can 
find for these systems an intimate connection between their injection properties and their 
dynamical features, such rationale could perhaps be adapted to more realistic systems. 
Such a procedure has been successfully applied in the study of fluctuations of current for 
conservative systems 10 1 . 

n 

In this paper, which follows a former one pj, we study a one-dimensional model of 
dissipative system, which has the advantage to allow for an exact description. This model 
consists of a chain of spins subject to an asymmetric T = Glauber dynamics, and is driven 
out of rest by a Poissonian flip of one spin (see next section for details): this is one of the 
rare examples where a nontrivial stationary dissipative state can be entirely described. In 
fact, we generalized the symmetric model studied in [jj by allowing for an asymmetry in the 
diffusion dynamics. This system is, in comparison with real dissipative systems, ridiculously 
simple, but one can hope that such examples would give ideas to interpret real experiments 
or to explain measurements on other variables, correlations, etc. . . For that purpose, this 
paper focuses much more on the physical content of the results than on the computational 
details, that are postponed in the appendix. More precisely, the observable we look at is 
the energy II provided to the system by the injection mechanism between t and t + r in the 
permanent regime. For large r, the probability distribution function (pdf) of II obeys the 
large deviation theorem and the probability distribution function is entirely governed by the 
large deviation function / (introduced below). The procedure of integrating the observable 
of interest over time has at least two advantages. First, one can hope that this effective low- 
frequency filter fades away "irrelevant" details of the dynamics and provides information on 
large-scale, hopefully more universal phenomena at work; this statement has been proved 
correct in some cases 



10 



llj . Secondly, the experiments are always constrained by a finite 



maximal frequency for the sampling of the time series: in practice the typical sampling time 
is much larger than the fastest relaxation times of the system under consideration. As a 
result, the pdfs experimentally measured are necessarily related to time integrated variables. 
The large deviation function is a good representation of these pdfs in the case where the 
sampling frequency is small with respect to the dynamics of the bulk. 

The paper is organized as follows: in the next section we define precisely the model; 
sections llHl and HVl are devoted to the physical results given by our computations. In section 
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IVl we show that the main characteristics of the large deviation function of the injected 
power are explained quite well using a simple phenomenological model, which treats the 
correlations between the boundary and the bulk in an effective way. The appendix (section 
IVIII) gives in detail all the steps of the computation, based on a free-fermion approach of 
the intermediate structure factor. 



II. THE MODEL 



We consider a ID system of N + 3 (N — > oo) classical spins on a line, labelled from —1 
to N + 1. The values of the extremal spins s_i and sn+i are fixed (this choice makes the 
description in terms of domain walls easier, as explained in the appendix). The zeroth spin 
So is the locus where energy is injected into the system: the flipping of sq is just a Poisson 
process with rate A, independent of the state of the other spins. The spins of the "bulk" , 
from Si to syv_i are updated according to an asymmetric T = Glauber dynamics. 

The asymmetric T = Glauber dynamics is defined as follows : given < p < 1, the 
probability for a spin Sj to flip between t and t + dt is 

dt[l - Sj([l -p]s j - 1 +ps j+1 )], (1) 

which is illustrated in figure HJ Note that if p > q (resp. <), the domain walls are locally 



2pdt 



2qdt 



2dt 



Odt 



FIG. 1: The rates of the asymmetric Glauber dynamics 



drifted to the left (resp. right); forp = q = 1/2 we recover the system studied in [1| (with 



the difference, that contrarily to 



1| the system is not duplicated on each side of s ; this 
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simplification yields simpler calculations and a physics a bit easier to analyze). Note that 
the case p < 1/2, where the domain walls easily invade the system, is probably the most 
relevant one for a comparison with experimental devices of turbulent convection. 

These dynamics are dissipative but a non trivial stationary state is nervertheless reached 
thanks to the Poisson process on sq which injects continuously energy into the system (the 
injected energy is positive on average; however, negative energy injections are also possible 
fluctuations due to the bulk dynamics). 



III. THE MEAN INJECTED POWER 



The mean value of the injected power (e) can easily be calculated. It is given by (e) = 
A[Prob(s = si) — Prob(s = — si)] = A(s s i)- To compute Uj = (s Sj) (we are interested 
here in the special case j = 1), we notice that the quantity Uj obeys a closed equation in 



the permanent regime (see HI] for details): 

-(A + l)Uj + pU j+1 + qUj-i = (2) 

(let us recall that q = 1 — p) with the boundary conditions Uq = 1 and Uoo = 0. The 
determination of Uj is simple : the polynomial pX 2 — X(X + 1) + q has a unique root r less 
than one, and therefore Uj = rK The mean injected power 

w = .A+wg+ira (3) 

2p 

is plotted in figure [21 One can see that it is an increasing function of A and a decreasing 
function of p. This last point can be easily understood, as for higher p the domain walls are 
more and more confined near the boundary, enhancing the probability of negative energy 
injection. On the contrary, for low values of p, the domain walls invade the system rather 
easily: for p < 1/2, they are drifted away from the site of injection. This is at the origin of 
a large positive value for the average injection energy. 



IV. THE LARGE DEVIATIONS OF THE INJECTED POWER 

The main result of our paper is the computation of /(e), the large deviation function 
of the injected energy. This is a central observable associated with the long time (or low 
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FIG. 2: Mean injected power as a function of A and p. 

frequency) properties of the fluctuations of the energy flux in stationary systems. Let us 
call II the energy injected into the system between t = and t — r; typically II scales like 
r for large r. The ldf /(r) is defined as 



f{e) = lim r- 1 ln[Prob(n/r = e)] 



(4) 



However it is simply defined, this quantity is difficult to compute or analyze theoretically, 
as it involves the knowledge of the complete dynamics of the system, and measures the 
temporal correlations which develop in a nontrivial way in the nonequilibrium stationary 
state. 

Usually, one computes first the ldf g{pt) associated with the generating function of II: 



{e aU ) 



~ e 

T— >0O 



rg(a) 



(5) 



More precisely g(a) = lim^oo r l log(e an ). Then, f(e) can be obtained numerically 
solving the inverse Legendre transform 



f(e) = min (g(a) - ae) 



(6) 



The details of the computation of g(a) are postponed in the Appendix. The formula for 
g(a) (equation (1551) ) is not easy to interpret physically. We are thus in a situation where 
the exact result does not really highlight the underlying physics, and in particular does not 
make the long-time properties of the injection process particularly transparent. In order 
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to clarify this, we will follow a very pragmatic way: first we will sketch the different ldfs 
corresponding to different values of the relevant parameters (A,p) and raise some questions 
associated to them. In the next section, we will see that some simple phenomenological 
models account very well for the observed behaviours (these models were neither discussed 

n 

nor even evoked in 

In figure [3] (a), we show various functions f{e) for different values of the parameters A 
and function of e/(e). f(e) is maximum for e = (e), which is a generic property 

of ldfs. Clearly, the curvature at the maximum is a major feature of these curves, and is 
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FIG. 3: (Color online) Large deviations functions /(e) as a function of s/{e) for various values of 
the parameters p and A. (a) curvature not rescaled, (b) curvature normalized 

strongly dependent on the parameters (X,p). Writing f(e) = — 1[— f"((e))(e) 2 ](e/ (e) — 1) 2 + 
o[ej (s) — l) 2 , we see that the relevant quantity associated to the curvature, once e has been 
rescaled by (e), is a = —f"((e))(e) 2 = g' (0) 2 / g" (0) . The curves rescaled by the curvature a 
are plotted in figure [3] (b), where it is seen that the curvature and the mean energy, though 
of primordial importance, are however not sufficient to characterize fully the ldf: there is 
no clear collapse of the curves. The dependence of a with respect to p and A is plotted in 
figure H] (a). Its behaviour is remarkably similar to that of (e) itself. Before explaining this 
point (see next section), we note that the mean value of the energy injected up to time r is 
(II) = r(e); besides, the (squared) relative fluctuations of this quantity is given for large r 
by [(n 2 )-(n) 2 ]/(n) 2 = l/(ra). The ratio of these two quantities [(U 2 ) - (U) 2 } / (U) = a /(e), 
called the Fano factor, is plotted in figure H] (b): it is comprised between 0.5 and 1.3 for all 
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FIG. 4: (a) Factor a = —f"((e))(e} 2 and (b) Fano factor a/(e), as a function of A and p. 

values of the parameters (A,p), which shows that the correlation between a and (e), though 
clearly demonstrated, is a bit loose. Thus, a sound question is to ask why a and (e) are 
correlated, and what are the factors which limit or modulate this correlation. These issues 
will be discussed in the next section. 

Let us go back to the rescaled ldf in figure [3] (b). One can notice that all curves display 
a noticeable counterclockwise tilt with respect to the parabola. As for the Fano factor, this 
tilt seems to be constant, with some minor relative differences. To quantify this tilt, one 
writes the Taylor expansion of /(e) /cr up to the third order like: 



A simple calculation gives x = fi ,/ "(0)fi ,/ (0)/fi , "(0) 2 . This parameter quantifies the tilt and 
can be a priori positive or negative. The variation of \ with (p, A) , plotted in figure [3, 
shows a rather complicated dependence of x with respect to the parameters (in particular 
an absolute minimum for p ~ 0.6 and A > 0.5), but with always x £ [0-6, 1]. Both the global 
trend and the finer details raise natural questions: why is the tilt is always positive ? What 
does it mean concerning the physics of the system ? Why is the dependence on p and A so 
complicated ? The next section provides a phenomenological model which give satisfactory 
answers to these issues . 



\{e/{e)-l) 2 + \{e/{e)-lf + o{e/{e)-lf 



(7) 
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FIG. 5: Parameter \ as a function of p and A. The left graph is a contour plot of the surface. 
V. DISCUSSION AND COMPARISON WITH A SIMPLIFIED MODEL 

In this section we compare the results obtained above to a very simple model, in order 
to see which global mechanisms are at work. 

We consider an oversimplified version of our system, called in the following "pure Pois- 
sonian model" or PPM. In this model, the injection is a Poissonian emission (with rate p) 
of domain walls (d.w.) into the system and the energy is incremented by one each time 
a d.w. is emitted. In this case, on has (II) = rp and [(II 2 ) — (n) 2 ]/(II) 2 = l/(rp), that 
is (e) and a are both equal to p. Thus, in our system, the global similarity between the 
two quantities, i.e. cr/{s) ~ 1 is not fortuitous, it is in fact a signature of the approximate 
Poissonian structure of the injection. 

Conversely, the violation of the relation a /(e) = 1, plotted in figure H] (b), is interesting, 
as it accounts directly for the coupling of the Poissonian injection and the structuration of 
the system near the boundary. In order to clarify this coupling, we can extend slightly the 
PPM to account for the variability of the Fano factor, by considering that in an "effective" 
energy injection, not only one domain wall is concerned, but in fact an average number 
of domain walls n& w . For instance, for p ~ 1, where the domain walls are confined at the 
boundary, the only way for the system to absorb energy is the following rare event: a domain 
wall is created between sq and si, it translates to the right (limiting factor), and then comes 
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back to annihilate with another entering domain wall. This "injection event" is so rare, 
that two such events are necessarily far apart from each other, and the statistics of these 
events is actually Poissonian. In fact, one can see that the effective number of domain walls 
associated to one event is two instead of one. Indeed, if one generalizes the PPM to emit 
ndw domain walls per event, one gets a / (e) = l/n^w' figure H] (b) gives Ud w — 2 in the p ~ 1 
region as expected. 

Another region is simple to analyze: for A ~ 0,p < 0.5 (the domain walls invade the bulk), 
the inner dynamics is so slow that the effective emission of domain walls invading the bulk 
is Poissonian; virtually no domain wall is reabsorbed by the boundary. One understands 
that this scenario breaks down rather abruptly when the drift is directed towards sq, which 
explains the singular behaviour of the Fano factor at (X,p) = (0,0.5). To summarize, the 
fact that cr/(e) < 1 for the most part of the parameter range illustrates the cooperative 
character of the energy injection. 

However, in the special case of large A, small p, one also expects a Poissonian behaviour, 
determined in this case by the natural time of the bulk dynamics: for very quick flipping of 
the spin sq, and p > 0, the emission of a domain wall into the system is only limited by the 
move of the first domain wall to the right. 

The PPM, even modified by the parameter n<j w is unable to account for the region where 
the Fano factor is larger than one (namely this regime of large A, small p): it is difficult to 
imagine an effective Poissonian emission of an average number of domain walls less than one. 
Moreover, if one also considers also the tilt parameter x, the disagreements are stronger, for 
it can be easily shown that the PPM (with n<j w allowed) yields x — 1> irrespective of the 
value of ridw- We conclude that if the global trend of a positive tilt is again a signature of the 
approximate Poissonian nature of the domain wall injection, the model is a bit too rough 
to account for the observed subtleties (except for regions where x — 1> which correspond to 
the cases commented above). 

In order to get a finer description of the phenomenology, we can add a new parameter in 
the PPM model. Instead of assuming a Poisson process for the emission of domain walls, we 
assume a Bernoulli process [ljj]: the time span [0, t] is divided into t/ At intervals of length 
At, during which n^w domain walls can be emitted with a probability pAt. One thus takes 
into account a possible waiting time after an emission event during which no other event is 
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on average allowed. For this model, one easily shows that 

(e) = pn dw (8) 

= 77^ — rr ( 9 ) 

' ^ 1 n dw (l - pAt) v 7 

1 - 2pAt . , 

* = T^ (10) 

We remark that now the Fano factor can reach values less than one. This is the case for 
small values of A and p < 1/2, where n dw is certainly one: here the Poissonian character 
of the process is imposed by A, but there can be a waiting period after a flipping of A, due 
to the finite time required for the bulk dynamics to remove the domain wall from its first 
position. 

We remark also that the x factor of the Bernoulli model is always less than one, exactly 
like in the real system. It confirms also our previous interpretation for the case A ~ and 
p < 1/2: a deep decrease of x is observed for increasing values of A, which is associated in 
the Bernoulli model with an increase of At. By the way, we can extract from the preceding 
equations the effective parameters n dw , At and p, knowing (e), a and x from our numerical 
computation: 

At = ^—^ (11) 
a 

n dw = (a/(e)r\2- X ) (12) 

P=— x (13) 

In figure El d [8] (a) and[8](b), we see the values of n dwi At, p and pAt respectively, extracted 
from the results for the spin model. The interpretation of figure [6] is obvious: as expected, 
the average number of domain walls nd w stays close to one for p > and all values of A, 
and reaches 2 for p — 1, the intermediate p values corresponding to a crossover region. 

The quantities p _1 and At are both related to the natural timescale of the effective 
injection process. The main difference between them is that p~ l is effectively the rate of 
the equivalent process, whereas At is somehow a "waiting time" during which two injection 
events have little chance to occur consecutively. At is plotted in figure [7J 

Figure El (a) shows as expected that the injection process is very inefficient for A ~ and 
also for p ~ 1; obviously, this curve is qualitatively related to (e), since the efficiency of the 
injection has immediate consequences on the mean injected power, but it is interesting to 
note that p is by no means constructed from (e), but from cumulants of higher order. 
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FIG. 6: Effective parameter function of p and A. 




FIG. 7: Effective parameter At as a function of p and A. 

Figured] (b) shows that the system is really Poissonian in the regions (A > 0,p < 1/2) 
and p < 1, despite the fact that At can be large (fig. [7]). Note also the vicinity of A = 
and p > 1/2, where something interesting happens: two factors that elsewhere favours the 
Poisson character of the process, namely A small and p > 1/2, are simultaneously at work 
here and act again each other. The results of this "collision" is that the process is clearly 
not Poisson for p around 0.8, due to huge values of At (see fig. [Tj); this is also a transitional 
region from a Poisson process with one domain wall to a Poisson process with two domain 
walls. 

Finally, when p < 1/2 and A is away from zero (this is the region where the Fano fetor 
is larger than 1), we find a non Poissonian process with At of the order of 10% to 15% 
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FIG. 8: Effective parameters (a) p and (b) pAt as a function of p and A. 

of p~ l . For instance, when p = 0, p increases as A increases: there is a crossover from a 
A-limited regime to a bulk-limited regime for which a waiting time is observed. In this case, 
the probability of two close injection events is weakened because an injection event uses two 
spin flips: one flips of s and then a flip of Si (for p = 0, s 1 can flip only if Si = — s ); thus 
the probability of two events within St goes like A 2 (<5t) 4 instead of (pSt) 2 for a pure Poisson 
process (PPM). This explains the emergence of the waiting time. 



VI. CONCLUSION 



In this paper, we have presented a one- dimensional model of a dissipative system, a half 
infinite chain of spins at T = in a Glauber dynamics with a drift toward or away the 
boundary, sustained in a nontrivial stationary state by an injection mechanism, namely 
the Poissonian flipping of the boundary spin. We computed exactly the large deviation 
function of the injected power and subsequently the first three cumulants of its probability 
distribution, which account for the mean value of the injected power, its fluctuations and 
the skewness of the fluctuations. Using a simple phenomenological model and its refined 
version, we have shown that it can account for the main physical characteristics of the 
injection process very convincingly, allowing for a relevant physical interpretation of the 
variations of the three cumulants with the parameters (X,p) (rate of Sq flipping, magnitude 
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of the drift), in terms of an effective rate of emission of energy "quanta", an average number 
of domain walls in each quantum, and a possible waiting time after an injection event. 

We can hope that this phenomenology could give an interesting scheme to interpret some 
experiments, where the same kind of injection mechanism is more or less reproduced. For 
instance, in a turbulent experiment, unpinning of vortices created near a moving boundary 
could be a process of energy injection suitable for the description framework that we propose 
here. We can also think of the bubble regime in the ebullition process, where the main part of 
the energy transfer occurs via unpinning of vapour bubbles. We hope that some experimental 
results 15j could find a simple interpretation in the kinetic description that we give here. 

Finally our mathematical calculations show that such simple out-of-equilibrium models 
are integrable: this opens the way to more generalizations. 



VII. APPENDIX: FERMIONIC APPROACH TO THE TIME-INTEGRATED 

INJECTED POWER 

It is useful to describe spin systems in the dual representation of domain walls : between 
the site j and j + 1 is located the possible domain wall labelled j for j = —1, ... AT. The 
state of the system is thus characterized by C = (n_ 1 , . . . ,n^), where the n« are either 
(no domain wall) or 1. There are 2 N+2 possible states in this representation; note that the 
domain wall n_i does not play any role, but is required to make the fermionic description 
tractable. The dynamical equation for the probability is given by 

N 

d t P(C) = X[P(Co) - P(C)] + 5>(C,MCj - C) - P(C)w(C -> C s )\ (14) 

i=i 

where Cj holds for the state C whose domain wall variables rij and Uj-\ have been changed 
(according to n — > 1 — n). The T = asymmetric Glauber dynamics corresponds to 

w(Cj — > C) — 2[1 — prij — qrij_i] (15) 
w{C — > Cj) = 2[prij + qrij-i] (16) 

where rij and rij_\ are the variables associated with the state C (we use this convention 
hereafter) , and q = 1 — p. 

We consider that each domain wall contributes as an excitation of energy 1 to the global 
energy of the system. We are interested in the energy II injected into the system up to time 
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t by the Poissonian injection. Following 13[, the route to this time integrated observable 
begins with the consideration of the joint probability P(C,H,t), the probability for the 
system to be in the state C at time t having received the energy II from the injection. The 
dynamical equation for this quantity is readily 

d t p(c, n) = \{p(c , n - i)n + p(c , n + i)(i - n ) - p(c, n)} 

N 

+ J2l p (Cj, K)w(C 3 - C) - P(C, U)w(C -> %)] (17) 



We define next the generating function of II as 

CO 

F(C) = eaUp ( c ^ u ) 



(18) 



n=-oo 



This quantity, summed up over the states, yields the generating function (exp(aTI)) from 
which one derives its ldf g(a): 



(e Qli ) ~ e 

t— >oo 



tg(a) 



(19) 



This ldf g(a) is closely related to f(p), the ldf of the probability density function of II, as 



they are Legendre transform of each other 
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Prob(II/t = p) oc exp(tf(p)) 



(20) 
(21) 



f(p) = min (g(a) - op) 

a 

Let us write the dynamical equation for F : 
d t F(C) = A [e a F(C )n + e" a F(Co)(l - n ) - F(C)] 

N 

+ 2 ^[F(C y )(l - pn j - qnj-i) ~ F(C)(pnj + qn^)\ (22) 
i=i 

The function g(a) can be expressed in terms of the linear operator acting on the "vector" 
[F(C)]c in the r.h.s of (T221 : it is in general its largest eigenvalue. 

Our problem belongs to the category of the "free-fermions" problems, for which a diag- 



onalization o 
described in 



the dynamics into independent "modes" can be achieved. The procedure is 
1], with references therein. In our case, the operator in the r.h.s of equation 
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( 122]) can be turned into the following fermionic operator : 
H = A[e a cLi4 + e"°c c_i + e a c\c. x + e- Q c t „ 1 c - 1] 

N 

+ 2 ^[cjCj_i + qc)c^i + pct_ x c,- - pc^cj - qc^cj-x] (23) 

3=1 

A symmetrisation procedure is a prerequisite to solve the problem. We define a priori the 
change of variables (note that the c remain fermionic variables) 



c_i = e a c_i , clx = e a cLi 



Vj ^ l,Cj = MjCj , c] 



■4 



where the Uj are real quantities to be defined. The choice 



u i 



3/2 



j>0 \ y p / 

leads to the symmetrized expression (we omit the tildes immediately) 



H = A[e 2Q cl lC J + e- 2Q c c_i + cjo-i + c T _ lCo - 1] 

JV 

+ 2j>^ 1 / 2 

3=1 

N 

= £ 



,,t. 



,t 



n,m=— 1 



A 



(24) 
(25) 



(26) 



"CjCj-i + yJ-pqc'jCj-x + y/pqc^Cj - pc]cj - gcj_ 1 c i _i] (27) 

(28) 

where A is a (AT + 2) x (AT + 2) tridiagonal, real and symmetric matrix, and B and -D 
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(N + 2) x (N + 2) antisymmetric real; they are defined by 

/O A 

A — 2q 2^/pq 



A = 



2y/pq —2p J 



( 



B 



Xe 2a 
-Xe 2a 



\ 







D 



( -Xe~ 2a 
Xe~ 2a 



-2i/3 



2v% 



-2v\ 







-2i/i 



-A 



-JV 



This Hamiltonian is diagonalizable, that is, it can be written 



(29) 



(30) 



(31) 



(32) 



where the £ g are fermionic operators linearly related to the Cj and the eigenvalues A q are 
the eigenvalues with a positive real part (we could have chosen the other half as well, see 
below) of the matrix 



M n 




(33) 



(The details of this procedure are exposed in [l|; note that the lack of translational invariance 
prevents the use of a Fourier transformation). 
The eigenvalues of H are thus given by 



N - A 



(34) 
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where the e q are ±1. In particular, the largest eigenvalue of H reads 

g(a) = ^Re(A q )-N-\ (35) 
<? 

= — £ d/i/i^4 -N-X (36) 
Am J Xo(/i) 

where Xo is the characteristic polynomial of M and the contour of integration is diverging 
half circle leant on the imaginary axis, with its curved part pointing toward the region 
Re(/x) > 0. 



A. The characteristic polynomial : introduction 

The problem is now equivalent to finding the characteristic polynomial of Mq. We can 
take advantage of the emptiness of B. We define E(ji) = (A + fi)~ 1 D(A — /f) -1 . Multiplying 
//Id - M by 

(/i-A)" 1 
(/i + A)~ 1 D(fi - AY 1 (n + A)- 1 

we see that 

xM = Xa{v)Xa{-h) det(l + BE{p)) (38) 

where Xa(^) = det(A — fi). Besides, 

det(l + BE) = (1 - Ae 2a £_i, )(l + Ae 2o £ ,-i) + A 2 e 4a £ ,o£-i,-i (39) 
= [1 + Ae 2Q E ,_i(-/i)][l + Ae 2Q £ ,-iGu)] + A 2 e 4Q £ , (/i)£_i,_i(/x) (40) 

where we exploited the fact that E T (fi) = —E(—fi). Note in passing that the symmetry 
of this expression with respect to fi — > —fi, is here explicit, as the Ejj are antisymmetric 
functions of \i. 




B. Some minors of //Id + A 

We term Aj, (j = 0, . . . , N + 1) the determinant of the minor of (// + A) obtained by 
keeping the (N + 1 — j) x (N + 1 — j) matrix located at the bottom right side of (/x + A) (one 
adopts the convention Ajv+i = 1). Note that A N = —2p + /i, A N _x = (/i — 2)(/x — 2p) — 4pg, 
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and that we have that det(/x + A) = /xA — A 2 A X . We have also an explicit formula, valid 
for j ^ 0: 

1 



A 



1 — Apqx 2 , 



(2qx + + l)^^^" 1 - 2qx+(2px+ + l)(Apqx+)- j+N+1 (41) 



where x + is conventionally the root of the polynomial — ApqX 2 + (/x — 2)X — 1 = with the 
largest modulus. Note that x + and Aj depends on \i. Later on, we will denote x_ = /x) 
and Aj = Aj(— /x); let us stress here that x + and x_ are roots of different polynomials. 

C. The characteristic polynomial : explicit calculation 

From the definition E = (A + fi)~ 1 D(A — it) -1 , one gets after some computations 



W 



1 N 

-5>^'[A i+1 A--AT +1 A,] (42) 



2q 
H 7=1 



^o,-i = JZrrT—^TT—A [^ 2a M^ ~ A 2 AD + 2„W] (43) 



det(A + /i) det(A 


- fJ-) 


-A 2 




det(A + /i) det(A 


- fi) 


-n 





= det(.4 + „)det(.4-,) [ (AlA » " A ' A ° )C "° + 2W T (44) 

^ det^de^-,) (45) 
whence one deduces 

Xo (ji) = -/i 2 A A + 2A 2 /x[A A x - A Ar] + 4A 2 /xe 2c W (46) 



Let us analyse W. Using equations (1421) and (JHj), one easily shows that PU has 5 different 
terms, respectively proportional to ((4pg) /4pqx_) N , (x-/4pqx + ) N , 

and (4g 2 ) Ar . From the definition of x + and x_, one has always 4pg|x-t| 2 > 1. This shows 
that the first term always dominates all but the last. A slight issue arises here, for the last 
is not always the least: for q < 1/2, the first is still dominating, but for q > 1/2 this is 
not the case for all values of /x in the right half plane V = {/x; Re(/x) ^ 0}. Let us term J 
the zone in V where 4g 2 ^ (4pq) 2 \x + x_\. Two key features of J are that (i) it is bounded 
(compact) (ii) it crosses the vertical line Re( /x) = only at one point, /x = 0. To prove that, 
we remark that on that line, x_ = x* + (complex conjugate); moreover x+(/x = 0) = l/2p 
and x+(—iy) = x + {iy)*: the maximum principle leads to the conclusion that \x+(iy)\ is a 
function of y, minimum at y = 0. 
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As a result, the contour of integration in equation (1361) can always be chosen such that, 
except for the single point \i — 0, it does not cross the region J (it encloses it anyway). In 
that case, the thermodynamic limit can be safely taken for all values of q, and leads to the 
complete vanishing of the term proportional to (4q 2 ) N in the result, dominated by the first 
one. This mathematical argument yields a great simplification, as one can consider that at 
the thermodynamic limit, W is always dominated by the term oc [{Apq) 2 x + xJ\ N and throw 
away the others. 

According to the preceding discussion, we are left with 

W= f A+A - (16p^x + x_) N (x_ - x + ) (47) 

A ± = 2gX±(2 ^ + 1) (48) 
Apqx± — 1 



(49) 



Similarly, we can write for j ^ 1 



A = A + (4pqx + )-i +N+1 (50) 

N— >oo 

As regards A , we have A = (— 2q + /u)Ai — 4pqA 2 . Thus, 

A = A + {4pqx + ) N [-2q + n-l/x + ] (51) 

N— >oo 

= A+(Apqx+) N x2p(l + 2qx+) (52) 

As a result, we get 

XM = A + A_{{Apqfx + xJ) N 

- 4pV(l + 2gx+)(l + 2qx.) + 8pq\ 2 fi{x- - x+) + 8p\ 2 fie 2a 



Ap 2 x + x^ — 1 
(53) 



A + A- ( (4pq) 2 x+x_ ) N Apu 

x ( - pf i(l + 2qx+)(l + 2qx.) + 2A 2 (x_ - x+) A fJ^~ + f + 2A 2 (x_ - x+) T ^— 



4p 2 x + x_ — 1 4p 2 x + x_ 

(54) 

This expression can be transformed in the following way: we can demonstrate the relations 

1/1 ^ (55, 



4pqx + X- — 1 2yU x + 

1 + 2qx + ){l + 2gx_)(4p — 1) = // x+x- (56) 

1 — Qpqx—X-^ 
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(for the first, multiply the lhs by (x + /x^ — 1) 1 ; for the second, use the fact that x± are 
roots of 2 nd de gree polynomials). Thus we can write 



Xo(fJL) = A + A^({Apq) 2 x + x_) N Ap 2 

( 2X 2 u e 2a - 

x (4A 2 - fi 2 ) (1 + 2qx + )(l + 2gx_) + -^(x_ - x+)—^ 

V V Ap z x + x 

= A + A_((Apq) 2 x + x_) N Ap 2 {AX 2 - fi 2 ) (1 + 2qx+)(l + 2gx_) 



AX 2 /p e 2a - 1 

x 1 + ,x a 1 I i (58) 



4A 2 — /i 2 4pgx + x_ + 1 

D. An integral formula for 5 

It is useful for the sequel to give explicit formulas for x + and x_ in the half plane 
V = {//; Re(/i) ^ 0}. A careful inspection shows that x + is given by 

( (Spg)- 1 X (/i - 2 - v/(^ - 2) 2 - 16pq) if Re(^) G [0, 2] 

= < ) < (59) 

[ (8m)- 1 x U - 2 + v/(/i - 2) 2 - 16pq) if Re(/x) > 2 

It must noted that, contrary to appearances, x + is analytic on the line Re(/x) = 2. It has 
however anyway a branch cut, localized on the segment \i G [2 — A^/pq, 2 + A^fpq]. 
The behaviour of x_ is entirely different: 

x_(//) = -(8m)- 1 x + 2 + V(^ + 2) 2 - Wp~q) (60) 

and is analytic on V. 

Let us go back to formula fl36l) . We see that only the logarithmic derivative of xo is 
involved, so we can handle the different terms of (1581) separately. For sake of clarity, we 
define 

rl ^ll^m (61) 

over a contour (followed counterclockwise) in V large enough to encircle all the singularities 
of/- 

• the terms x±: as x_ it is an analytic function of \i over V, we get I[x^} = (obviously 
neither x_ nor x + can go to zero). The term I[x+] gives a nonzero contribution 
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due to the branch cut of x + . We remark that on it, \x. 
counterclockwise the circle of radius 2^/pq. Thus, 



I[x 



N I dZ <* 1 OX 

— & —{Apqz + Z - 1 + 2) 

J\z\=2y/pq Z 



2^/pq and x + describes 



(62) 



N 



The term (4A 2 - /i 2 ) yields J[4A 2 - /i 2 ] = A. 

We show easily that A±(l + 2qx±) = ^jix\/ (Apqx\ — 1). We conclude easily that 
ILA_(1 + 2qxJ)} = 0, for 4pgx 2 — 1 never vanishes. As regards IL4.+(1 + 2qx + )}, a 
transformation similar to (|62[) gives also 7L4 + (1 + 2qx + )] = 0. 



Finally, from these results, we see that all terms but the last in (1581) cancel with constant 
terms in fl36l) . To give the final result a convenient form we remark that the contour of 
integration can be make infinite, that the semicircular part gives a vanishing contribution 
to the result, and that on the vertical line Re(/i) = 0, x + = x*_. We can thus write 



9{oc) = - 

71 



dy log 1 



A 2 /4p 



,2a 



A 2 /4 + y 2 ip(y) + 1 



tp(y) = 4pq\x_(Uy)\ 2 = (4pq) 1 x 2iy + 1 + y/(2iy + l) 2 - Apq 



(63) 



(64) 



We verify immediately that g(0) = as expected. We can also check that </(0) = (n) 
i2 r °° dy 



9 '(0) = ™ I 
pn J_ 



-A 



oo 

2 roo 



(A 2 + 4y 2 )(^(y) + r 
dy ( 1 



dy 



1 



2pn J„ oc A 2 + Ay 2 \x + (4iy) x~(4iy) 

\2 poo 

--Re 

pn 



A 



— (A + l- v / (A + l) 2 -4 M ) 



(65) 
(66) 
(67) 



, A 2 + Ay 2 x_(4iy) 2p 

We mention (without computations) also the result we would obtain, if we had considered, 
like in two half lines of spins connected to s (i.e. spins numbered s_i, s_2, . • •)• In this 
case, despite the fact that the two subsystems are connected only via the Poisson spin so, 
they are nontrivially coupled to each other, and the corresponding large deviation function 
of the cumulants reads 



A 2 /4p 



3 2a - 1 

9 {a) = - I dy\og I 1 , ' ' ' ~ 

n j V A / 4 + y v(y) + 1 



,2a 



2 + 



- 1 



p(ip(y) 



(68) 
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We see clearly that there is no simple correspondence between the half line model and the 
two half lines model: g is multiplied inside the logarithm by an a-dependent term. 
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